FIRST ORDER CONTINUOUS MODELS OF OPINION FORMATION* 

GIACOMO ALETTI+, GIOVANNI NALDI* , AND GIUSEPPE TOSCANI§ 

Abstract. We study certain nonlinear continuous models of opinion formation derived from 
a kinetic description involving exchange of opinion between individual agents. These models imply 
that the only possible final opinions are the extremal ones, and are similar to models of pure drift 
in magnetization. Both analytical and numerical methods allow to recover the final distribution of 
opinion between the two extremal ones. 
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1. Introduction. This paper is devoted to the analysis and large-time behavior 
of solutions of the equation 

§=7^((l-z 2 )(*-™W)/) (1.1) 

where the unknown f(x, t) is a time-dependent probability density which may rep- 
resent the density of opinion in a community of agents. This opinion varies between 
the two extremal opinions represented by ±1, so that x el = [—1,1]. The constant 
7 is linked to the spreading (7 = —1) or to the concentration (7 = +1) of opinions. 
In (1.1) m(t) represents the mean value of /(■,£), 

m(t) = jj{x,t)dx, (1.2) 

and its presence introduces a nonlinear effect into its evolution. The value of the 
constant 7 is responsible of completely different effects. In case 7 = — 1, this equa- 
tion, from now on called nonlinear decision equation, arises in the study of opinion 
formation, and has been introduced in [16], in connection with the quasi-invariant 
opinion limit of a model Boltzmann equation for the kinetic description of opinion 
formation involving exchange of opinion between individual agents. The equilibrium 
state is given by two Dirac masses located at the extremal points of the interval I. In 
case 7 = +1, the equation represents a simple type of nonlinear friction equation, and 
its effect is to concentrate the solution. In this case, the equilibrium state is given by 
a Dirac mass located at some point of the interval T. 

When 7 = — 1, related equations have been introduced recently. In [11] a linear 
equation describing the pure drift in magnetization has been derived as the mean field 
limit of the Sznajd model [14] in case of two opinions. There the (linear) first-order 
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partial differential equation coincides with equation (1.1), where it is assumed that 
m(t) is identically zero. This simplification allows for an analytical treatment, since 
it is possible to obtain the exact solution and control the rate of decay towards the 
equilibrium for all values of 7. If 7 = +1 related models of one-dimensional nonlinear 
friction equations have been considered in the study of granular flows in [9, 15], in 
connection with the quasi-elastic limit of a model Boltzmann equation for rigid spheres 
with dissipative collisions and variable coefficient of restitution. 

Microscopic models of both social and political phenomena describing collective 
behaviors and self-organization in a society have been recently introduced and ana- 
lyzed by several authors [5, 8, 10, 12, 13, 14, 19, 20]. The leading idea is that collective 
behaviors of a society composed by a sufficiently large number of individuals (agents) 
can be hopefully described using the laws of statistical mechanics as it happens in 
a physical system composed of many interacting particles. The details of the social 
interactions between agents then characterize the emerging statistical phenomena. 

Among others, the modeling of opinion formation attracted the interest of a in- 
creasing number of researchers (cfr. [4, 8, 10, 14] and the references therein). Mean 
fields model equations have been proposed in [1, 11]. These equations are in gen- 
eral partial differential equations of diffusive type, that can in some case be treated 
analytically to give explicit steady states. A kinetic description based on two-body 
interactions involving both compromise and diffusion properties in exchanges between 
individuals has been proposed in [16]. Compromise and diffusion were quantified by 
two parameters, which are mainly responsible of the behavior of the model, and allow 
for a rigorous asymptotic analysis. In a compromise dominated regime, the resulting 
equation is our equation (1.1). 

The mathematical methods we use are close to those used in the context of kinetic 
theory of nonlinear friction equations [7] , and made popular by the mass transporta- 
tion community [18]. 

The paper is organized as follows. In the next section we introduce the main 
properties of the model, which justify the treatment in terms of a suitable weak for- 
mulation. The qualitative analysis is given in Section 3. The large-time behavior is 
considered in Sections 4-6. It is shown that the problem can be solved in sufficiently 
high generality only in the case of concentration (7 = 1). This lack of generality in 
the analytical treatment of the large-time behavior of the solution in the spreading 
of opinion justifies the numerical treatment of the equation. The numerical approxi- 
mation is included in Section 7. 

2. Main properties and weak description. As briefly described in the in- 
troduction, equation (1.1) describes the evolution of a probability density which rep- 
resents the density of opinions in a community. For all values of the constant 7, we 
will show that the time-evolution driven by this equations leads the density towards 
a equilibrium state that is described in terms of two Dirac masses (7 = —1) or to a 
unique Dirac mass (7 = 1). Having in mind that the equilibrium solution to equation 
(1.1) is given by Dirac masses, any convergence result towards equilibrium holds in 
weak*-measure sense. The recent analysis of [7] of the nonlinear friction equation 
introduced by McNamara and Young [9], suggests that a suitable way of treating 
equation (1.1) is based on a rewriting of this equation in terms of pseudo-inverse 
functions. It is immediate to show that the drift operator on the right-hand side of 
(1.1) preserves positivity and mass, 




(2.1) 
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Then, given a initial datum which is a probability density (nonnegative and with unit 
mass), the solution remains a probability density at any subsequent time. Let F(x) 
denote the probability distribution induced by the density f(x), 

F(x)= [ f(y)dy (2.2) 

J ( — 00, x] 

and let p denote the distribution on R associated to F. Since F(-) is not decreasing, 
we can define its pseudo inverse function (also called quantile function) by setting, for 

pe(o,i), 

X»(p) = X F (p) = M{x : F(x) > p} 

Equation (1.1) for f(x, t) takes a simple form if written in terms of its pseudo inverse 
X(p, t). Theorem 3.1 shows in fact that the evolution equation for X(p, t) reads 

= _ 7 {X (p, t) m{t)) (1 - X 2 (p, *)), (2.3) 

where now p € (0, 1). Note that if we assume F to be absolutely continuous with 
respect to x and strictly increasing, then Theorem 3.1 reduces to elementary compu- 
tations. In (2.3) 



m(t)= / X(p,t)dp. (2.4) 
Jo 

Let us set 7 = — 1 (spreading). Then, the weak form (2.3) clarifies the evolution of 
X(p, t) and the role of m(t). In fact, if X(p,t) > m(t), X(p,t) increases towards 1, 
while X(p,t) < m(t) implies that X(p 1 t) decreases towards —1. Hence, the mean 
opinion m(t) represent a barrier for the density of opinions to move towards one of 
the two extremal opinions. The fact that the mean opinion varies with time makes 
the nonlinear problem harder to handle with respect to the linear problem considered 
in [11] where the barrier is fixed equal to zero. 

Among the metrics which can be defined on the space of probability measures, 
which metricize the weak convergence of measures [21], one can consider the LP- 
distancc (1 < p < 00) of the pseudo inverse functions 

d p (X,Y)= (j\x(p)-Y(p)\Pd p y /P . (2.5) 

In what follows, we'll use the usual identifications 

d p (X,Y) = dpifxjy) = d p (F x ,F Y ) = d p {p x ^Y) , 

where px {^y), Fx (Fy) and fx (fy) denote the distribution, the cumulative function 
and the density associated to X (Y), respectively. By this identification, as one can 
see [7, 17, 18], d2(F,G) is nothing but the Wasserstein metric [17]. 

In addition to nonlinear friction equations arising in the modelling of granular 
gases [7], the strategy of passing to pseudoinverse functions has been recently ap- 
plied to nonlinear diffusion equations of porous medium type [3], and to degenerate 
convection-diffusion equations [2] . This rewriting of nonlinear diffusion equations has 
been shown to be useful in order to obtain simple explicit numerical schemes that 
satisfy a contraction property with respect to the Wasserstein metric [6] . 
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3. Existence, uniqueness and well position of the problem. In this sec- 
tion we will study the initial value problem for equation (1.1), with initial density 

f(x,t = 0) = f (x), xel. (3.1) 

As before, we will denote by X (p) the quantilc function corresponding to /o, so that 

X(p,t = 0)=X o (p)=mf{x:F o (x)>p}, pe[0,l]. (3.2) 

The equivalence between equations (1.1) and (2.3) is contained into the following. 

Theorem 3.1. There exists a weak solution of (1.1)— (3.1) if and only if there 
exists a solution of (2.3) -(3.2). 

Proof. Suppose first that there exists f(x,t) which solves (1.1)— (3.1). Then m(t) 
is a differentiable function of time. Let y{t) be the maximal C 1 solution of the Abel 
differential equation: 

y' = -7(1 - v 2 )(y - m (t)) , 33) 

2/(0) = Vo 

where, for any y € [—1,1] we denoted by yo a C 1 -extension of yo to K. We have, in 
weak sense, 



d_ 

dl 



f(x,t)dx= [ \^-( l { _ 00ty{t)] (x))f(x,t)+ l { _^ v(M {x)^-f{x,t) 
(-oo.y(t)] Jvi l ot\ j at 

-j 

Jr 



dx 



(y\t)d yit) (x)+ 1 [S y{t) (x)(l-x 2 )(x-m(t))])f(x 1 t)dx 
r (3.4) 

= 0. 

Since X(p, t) < x p < f(y, t) dy, the first part of the proof has been shown. 

Now, let X(p,t) be a solution of (2.3)-(3.2). As a consequence of the properties 
of the solution to Abel's equation (3.3), given any initial datum X(p, 0) satisfying 

• X(p,0)e [-1,1]; 

• X(p, 0) is nondecreasing; 

• X(p, 0) is left-continuous. 

the same properties are preserved at any subsequent time t > 0. Hence, for any 
t, {X{p 1 t), p £ (0,1)} is the quantile function of a unique probability measure on 
[— 1, 1]. We have only to prove that (1.1) holds. This is a consequence of the change 
of variables formula. In fact, if h is a test function 



/ M x )-57/0M)c^ = -57 [ h(x)f(x,t)dx 
Jr at at Jr 

= l- t j\ { x p{ t))dp 



1 d_ 

, dt 



h{X p (t))dp 



= J o ti(X p (t))( - 7 (X(p, t) - m(t)) (1 - X 2 (p, t))) dp 

= I h'(x)(-f(x-m(t))(l-x 2 ))f(x,t)dx 
Jr 

= J h(x)^ (7 (x - m(t)) (1 - x 2 )f(x, *)) dx . 
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□ 

We call (JC,p) the (compact) set of probability distributions on [—1, 1] equipped 
with the p-Wasserstein distance. Note that all the p-Wasserstein distances on (JC, d) 
are equivalent. In fact, if q > p > 1, 

H-X"" 1 -X^\\ p < \\X^ -X"*\\ q < 2 1 - p / q \\X f * 1 -X^\\P/ q . (3.5) 

We will refer to JC as the topological space of probability distributions on [—1, 1] 
induced by any of this metric: the weak*-topology. Before searching for a continuous 
solution of (1.1) in JC, we state the following trivial lemma. 

Lemma 3.2 (Solution Abel). Let <j>(x,y) = -7(1 - x 2 )(x — y). Then 

|<H#i, 2/1) - 0(^2,2/2)1 < 4|xi - x 2 \ + \yi - y 2 \ ■ 
Moreover, if f is a solution of f'{t) = <p(f(t),g(t)) with snp\g(t)\ < 1 and f(0) G 
[-1.1], 

|/(*)-/(t)|<2|s-t|. 

We call solution of (1.1)- (3.1) any function p, G C°(R,JC) s.t. (1.1)-(3.1) holds. We 
have the following 

Theorem 3.3. For any probability density fo{x) in (3.1), there exists a unique 
function p, G C°(R, JC) such that, if f(x,t) denotes the weak derivative of the probabil- 
ity distribution jJ,(t), f(x,t) satisfies (1.1) with initial value (3.1). Moreover, for any 
t e R, the solution depends in a continuous way on the initial datum: the problem 
(1.1)-(3.1) is well-posed in C°(R,JC). 

Proof. [Existence] We prove the existence of a solution of the equivalent problem 
(2.3)-(3.2) (see Theorem 3.1) in a constructive way. More precisely, 
A we construct a sequence {X n ,n € N} which approximates a target solution; 
B by compactness arguments, we find a convergent subsequence X ni — ► X; 
C the limit X satisfies (2.3)-(3.2). 

Let [— T, T] be fixed. For any n G N, we subdivide [— T, T] into disjoint intervals 
of length R/2 N . Then we proceed as follows: 
Al we compute = X a (p) dp; 

A2 we solve (2.3) on [-T/2", T/2 n ] with mS n \t) = m , finding X^(p,t), t G 

[-T/2",T/2™]; 
A3 for any k = l,...,2 n - 1, 

• we compute 

m£2 = f X (n \p,±kT/2 n )d P ; 
Jo 

• we solve (2.3) on (kT/2 n , (k + l)T/2 n ] with m^(t) = m ( k n) and initial 
dataX(")(p,fcT/2 n ), finding X^ n \p,t), t G (fcT/2™, (k + l)T/2™]; 

• we solve (2.3) on [-{k+l)T/2 n , -kT/2 n ) with mS n \t) = m M k and initial 
dataX(")(p,-fcT/2"), finding X^ n \p,t), t G [-(k + l)T/2 n ,-kT/2 n ). 

We call /x^™) : [— T, T] — > /C the sequence of function with value in /C associated to 
X (n \ 

B For any n G N, it is possible to prove (by induction on k) that for any t G [— T, T] 
X^ n \p,t) G [-1,1] and m(")(i) G [-1,1]. As a consequence of Lemma 3.2, we have 



(3.6) 
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i.e., for any p e (0,1), {X^(p,-) : [—T,T] — > [— 1, l]} n eN is a uniformly equicon- 
tinuous sequence. A diagonal argument together with Ascoli-Arzela Theorem ensure 
the existence of a subsequence ni s.t. {X^ ni \p, •) : \—T,T] — ► [— 1, 1]}/ £ n converges 
uniformly on [-T, T] for each p <= (0, 1) n Q. 
Now, (3.6) implies 



/ \X^(p,s)-X^(p,t)\dp<2\t-s\ 
Jo 



(3.7) 



i.e., is a equicontinuous sequence with respect to the distance di on K. Then, 
Ascoli-Arzela Theorem again ensures the existence of a subsection of n ; (we call it n ; 
again) s.t. 

sup di(/i( n *>(t),/i( n, >(i)) < M{kM) — ► (3.8) 
sup |X( n *)(p,t)-X( n ')(p,t)| < N p (kAl) — » 0, VpeQn(0,l) (3.9) 

t£[-T,T] kAl^oo 

Now, let p € (0, 1) n Q be fixed. {X(" ! )(p, -)}zeN is a uniform convergent sequence 
of derivable functions converging to X(p,t) = limi X^ n '\p,t)). Left-continuity and 
monotonicity of {X(p, t),p G (0, l)nQ} extend the definition of X(p, t) to all p € (0, 1). 
C What remains to prove is 

• lim;m^)(i) = X(p, t) dp =: m(t); 

• X(p,t) is differentiable, and (2.3) holds. 

Note that, from (3.8), it follows immediately that \ ^ X {nk \p,t) dp- ^ X {ni \ Pl t) dp\ 
< M(k A h) and hence 

f X(p,t)dp- [ X < - n '1(p,t)dp <M(h). (3.10) 
Jo Jo 

By definition of m^ n \ 

m^(t) = m (n \[[2 n t]}/2 n ) = f X^ n \p, [[2 n i]]/2 n ) dp 

Jo 

where [[•]] is the integer part of • closer to 0. Therefore, 



in 



(*)- / X^(p,t)dp 



< 2 



-n+l 



(3.11) 



and hence m^ ni \t) is a Cauchy sequence on [—1,1]. Thus, there exists m(t) = 
\\m h m {ni \t). By (3.10) and (3.11) it follows that fh(t) = m(t) since 

»(t)- / X(p,t)dp < \m{t) - m^ ni \t)\ +2- ni+1 +M(h). 
Jo 

Now, let p e (0, 1) ("1 Q be fixed. For simplicity of notation, define H(n,p,t) := 

d 

—X( n \p,s) . Moreover, we define 
as 



H(p, t) := lim H(m,p, t) = - 7 (1 - X(p, t) 2 )(X(p, t) - m(t)). 
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The uniform convergence theorem states that —X(p, s) 

as 



= H(p,t) i{{H(n u p,t)} l& 



is a uniform converging sequence on (— T, T). To prove this, let k > I. Triangular 
inequality 

\H(n k ,p,t)-H(ni,p,t)\ < \H(n k , p, t) - H(n k , p, [[2 n, t]]/2 n, )| 

+ \H(n k ,p, [[2"'t]]/2"«) - H(n u p, [[2 n 't]]/2 n ')\ 
+ \H(n l ,p,[[2 n H]]/2 n >)-H(n l ,p,t)\ 
= A p (k, I, t) + B p (l, k, t) + A p (l, I, t) 

shows that we may prove that sup te r_ TT1 A p (l, k, t) + BM, k, t) + A p (l, I, t) — ► 0. 

As a consequence of (3.7) and (3.11), we have \mS n \s) - mS n \t)\ < 2~ n + 2 + 2\t - s\, 
which implies (see Lemma 3.2 and (3.6)) 

|ff (n, p, s) - H(n, p, t)\ = | (1 - JTW (p, S f)(X^(p, s) - m W («)) 

- (1 - X< n \p,t) 2 )(x( n \p,t) - m W(i)) 
< 4|X(")(p, s) - X (n) (p,t)| + |to (,i) (s) - m W(t)| 



< 10|i ■ 



-n+2 



and hence < 10-2-™ ! +2-™ fc + 2 . Now, let fe > I and / e {-2™'+l, . . . , 2"<-l}. 

Again, Lemma 3.2, (3.8) and (3.9) imply 

\H(n k , p, l/2 n >) - H(n u p, Z/2"')| < 4JV„(0 + M(J) 

and hence B p (k,l,t) < 4N p (l) + M(l). This completes the proof for p e (0, 1) n Q. 
Now, fixed j/o € [ — 1, 1] , let y(t) be the maximal C 1 solution of the Abel differential 
equation (3.3). Since X(p,t) > y{t) <=^- X(p,0) > 2/(0), Vp e (0,1) n Q, left- 
continuity and monotonicity of {X p (t),p € (0, 1)} extend the proof tope (0, 1). 

[Uniqueness and well-position] Let Y(p,t),X(p,t) be two solutions of (1.1)— (1.2). 
Denote by mx{t) and my it) the mean values of X and Y respectively at time t. 
Since |my(i) — mxit)\ < d\(X,Y) we have (by Lemma 3.2 and (3.5)) 

f t d 2 (p X (t),pyit)) = jj^ (Yip,t)-Xip,t)) 2 dp 

< 2^ \Yip,t)-X(p,t)\(4\Yip,t)-Xip,t)\ + \m Y (t)-m x (t)\) dp 

< 10 ■d 2 (px(t),pY it)). 

Gronwall's lemma completes the proof. □ 

4. Large— time behavior of solutions. Thanks to the uniqueness of solutions 
of Abel equation, we obtain the following 

Lemma 4.1. For any (^sel and p ^ p' e (0, 1), we have 

X(p,t) = Xip',t) X(p,s)=X(p',s), 

i.e. (1.1) does not create or destroy delta masses infinite time. 
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A direct consequence of the previous lemma is that the initial masses in +1, — 1 
and in (—1, 1) remain unchanged in time. Let us call themp +1 , p_i and 1 — (p+i+P-i), 
respectively. 

An important argument linked to the large-time behavior of solutions to nonlinear 
equations is the study both of conservation laws and of Lyapunov functionals. In 
addition to mass conservation, a second conserved quantity (when defined) is furnished 

by 

In addition to the conservation of both mass and T(t), equation (1.1) possesses a 
Lyapunov functional, simply given by the variance of the solution 

V(t) := J (X(p, t) f dp-(J o X(p, t) dp) 2 . (4.2) 

We give below an easy to check condition which insure both the the boundedness and 
the conservation in time of the functional (4.1). 

Lemma 4.2. Let log ((1 + X(p, 0))/(l - X(p, 0))) e L^O, 1). Then, for allteR 



is well-defined. Moreover, T(t) is differ entiable andT'{t) = 0. 

Proof. Since log ((1 + X(p, 0))/(l - X(p, 0))) e L^O, 1), then X(p, 0) e (-1, 1), 
Vp e (0, 1). Now (1 - x 2 ){x - 1) < (1 - x 2 ){x - m{t)) < (1 - x 2 )(x - 1), and hence 
the uniqueness of solutions of Abel equations imply 

X(p,t) e (-1,1), Vpe (o,i),VteR. 

Let G(p,t) = log((l + X(p,i))/(l -X(p,t))). Since \G t (p,t)\ < 2, we have 
. 3G t , Vp e (—1, 1),V£ e [— T,T]; 

• |G(p,t)| < |G(p,0)| + 2T and hence T(t) exists; 

• |G t (p,t)| < 2 and 2 e ^(0,1); 

and hence it is possible to differentiate under the integral sign, obtaining T'(t) = 0. □ 
The following lemma shows that the variance is a Lyapunov functional for equa- 
tion (1.1). 

Lemma 4.3. The variance V(t) of p(t) is a monotone differentiate function with 
values in [0, 1] . 

Proof. V(t) is clearly differentiablc, and 

V\t) = jj\x t (p)-m(t)) 2 dp = -2 7 j\l-X t (pf)(X t (p)-m(t)) 2 dp. 

In case 7 = — 1, V(t) is monotonically increasing, while bounded from above by 1. In 
fact, the maximum value of V is attained for p = (Si + S-i)/2. If on the contrary 
7=1, V(t) is monotonically decreasing, while bounded from below by 0. In this 
second case, the minimum value of V is attained for p — S a , a G [—1,1]. □ 

The most difficult problem linked to equation (1.1) is the study of the evolution 
of the mean m(t), and to the exact evaluation of its limit value rh = lim i ^ 00 m(t). 
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The knowledge of m is of primary importance, since in consequence of the structure 
of the limit state of the solution to equation (2.3), this value is enough to characterize 
completely the steady state. 

Remark. The difficulty comes out from the evolution of the mean m(i), which is 
given by the "non closed" equation 



In what follows, we make use of the previous results to extract information on 
the behavior of the mean m(t). 

Lemma 4.4. For any t, m(t) G [-^1 - V(t), y/l - V(t)]. 

Proof. Since V(t) = Jx 2 f(x,t)dx-(m{t)) 2 , then V(t) < 1 - (m(t)) 2 . □ 

Lemma 4.5. The function m! : R — > [-1,1] belongs to L 2 (R). Moreover, 
Blim^oo m'(t) = 0. 

Proof. It is sufficient to note that 



and hence m! G i 2 (K) by Lemma 4.3. Moreover, since ml is a Lipschitz function (in 
fact it is diffcrcntiable and m" < 10, cfr. (4.3)), it follows that \im t -,oo m'(t) = 0. □ 

5. Spreading of opinions. In this section we study the large behavior of (1.1) 
when 7 = — 1, leaving the opposite case 7 = 1 to the following section. 

Remark. If X(p, 0) = —X(l — p,0) (i.e., the initial distribution is symmetric), then 
from (4.3) follows that m(t) = for any subsequent time t > 0. If, in addition, 
X(pi,0) = X(pi +5,0) = 0, then X(pi,t) = X( Pl + S,t) = for all i > (i.e., any 
initial mass in is not moved away in time if the initial distribution is symmetric). 
In order to avoid these situations we will allow delta masses only in ±1: X(pi,0) = 



X(p 2 ,0) Pl = p 2 or (X( Pl ,0)) 2 = 1. 

Theorem 5.1. Assume X(p 1 ,0) = X(p 2 ,0) <=^> pi = p 2 or (X(p 1 ,0)) 2 = 1. 
Then, the limit distribution exists and it is given by two masses located in — 1 and 
+1. 

Proof. By Lemma 4.3 and Lemma 4.4, m(t) G [-^1 - V(0), y/1 - V(0)}, for all 
t > 0. Thus, if (X(p,t ) 2 > 1 - V(0), then (X(p,t) 2 > 1 - V(0), for all t > t . Since 

X(p, t) < -^/1-V(0) => ^X(p, t) < and X(p, t) > y/l - V(0) => ^X(p, t) > 
by Lemma 4.3 and Lemma 4.4, the two functions 



are monotone. We call p± ± — lim t _, 00 p£ 1 (t). (2.3) and monotonicity of Abel's 
solutions allows to state that 




(4.3) 





VpG [0,p\), lim X(p,t) = -1 
Vpe^i.l], }imX(p,t) = l. 



(5.1) 



(5.2) 
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Hence the limit distribution has two masses in —1 and +1. It remains to characterize 
what happens for the remaining p^_ 1 — p h _ x mass. Let us recall that 



x( P ,t)e [ - v/i - 7(0), y/i - v(0)] , Vpe^!,^), (5.3) 

By Lemma 4.5, there exists T > such that |m'(t)| < (F(0)/2) 2 for all t > T. 
By contradiction, suppose that there exists t > T and p g (f>-i)Pi) such that 
|X(p,i )-m(i )| > F(0)/4. 

Since |— X(p,i )| > |to'(*o)|, it follows that |X(p,i) - m(t)\ > V(0)/4, for all 



t > t a . Thus, (5.3) shows the contradiction: 

• ay- 

• X(p, t) is bounded. 
Therefore, 



j^X(p,t) is continuous; 
\i t x{ P ,t)\ > {V{Q)/2f if t > to; 



|X(p,t)-m(t)| < for all pe^), t>T. (5.4) 

Now, let F(x, y) = (1 — x 2 )(x — y) as in Lemma 3.2. Since F is differentiable, when 
x\ > x 2 , Lagrange theorem states that we can find £ g (xi, x 2 ) such that 



d 

F(x 1 ,y) - F(x 2 ,y) = {x 1 - x 2 )—F(x,y) 
Now, if 1 - x 2 > V(0) and \x - y\ < V(0)/4, we have 



^F(z, j/) - 1 - 3x 2 + 2xy > V(0) + 2x(y - x) > 



that is 



F( Xl ,y) - F(x 2 , y) > (an - ^2)^ , x ? < 1 - V(0) and | Xi _ „| < . (5.5) 

Let p^ < p 2 < Pi < P+i- Then, both (5.3) and (5.4) are satisfied for all t > T, (5.5) 
holds, and 

^- t (x( Pl ,t)-X(p 2 ,t)) > (X( Pl ,t)-X(p 2 ,t))^-, Vt>T. 

Since the two solutions are bounded, the only possibility is that X(pi,t) = X(p 2 ,t), 
for all £ > T , and for all (pi,p 2 ): p h L 1 < p 2 < p\ < p\ 11 which implies p h _ x = p\ by 
Lemma 4.1 and hypothesis. □ 

Remark. Theorem 5.1 may be read in terms of weak*-measure convergence 
f(x,t) - pfl 1 5- 1 (x) + (l-p?L 1 )6 1 (x). 

t—>oo 

In particular, since the support is compact, all the moments exists and will converge. 
We have the following 

Corollary 5.2. Assume X(p 1 ,0) — X(p 2 ,0) ^=^> Pl = p 2 or (X(p 1 , 0)) 2 = 1. 
Then there exists lim t ^ +tx) m{t) — rrioo = 1 — 2p^ 1 . 
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6. Concentration of opinions. Let us recall that the stochastic partial order is 
naturally given on /C. Let F(x), G(x) denote two probability distributions and Xp, Xq 
their pseudo-inverse functions, respectively. We say that F < G if F(x) > G(x),V(x) 
or, equivalently, if Xp(p) < Xg (p) , Vp € (0, 1). 

Lemma 6.1. The operator 

<f>(X) = -(X - m(X))(l - X 2 ) , 

is a monotone operator with respect the stochastic ordering. 

Proof. Assume X\(p, s) < X 2 (p, s), for all p e (0, 1). Then mi(s) < m,2(s) (they 
are equal iff the distributions coincide). Let p € (0, 1) be fixed. If X\(p, s) = X 2 (p, s), 
then X[(p, s) < X 2 (p, s). The continuity of X' is sufficient for the remaining part of 
the proof. □ 

Lemma 6.2. Let Xq in (3.2) be given. Then, there exists lim t ^ +00 m(t) — mac,. 

Proof. Let [a, b] be the class limit of m(t). Suppose a = — 1, i.e. liminf t m(t) = — 1. 
Markov inequality then shows that the limit distribution is a mass concentrated in 
— 1, and hence b = — 1. Otherwise, we may assume that m(t) e [—1 + S, 1 — 5], for all 
t > to and let po the mass not concentrated in ±1 at each time (recall Lemma 4.1). 
For all e > 0, po — e mass is in [—1 + e, 1 — e] at t = to- Therefore for all p € (0, 1) 
X(p, to) € [— 1 + e, 1 — e], and X(p, t) decades exponentially to m(t) with rate not less 
than (mm(6, e)) 2 . The large behavior of this process shows three delta masses: the 
initials in ±1 and the remaining in m(t). Stationary arguments imply the existence 
of moo. □ 

The steady state of the process can now be defined by the following 

Theorem 6.3. If (1 — Pi)(l — P-i) < 1 (i.e., if there are masses in ±1 at time 

t = 0) thenm^ =p 1 -p_ 1 . Otherwise, if log ((1 +X(p, 0))/(l - X(p, 0))) e £ x (0, 1) 

then 

_ cx P {T(0)}-l 
moc -cxp{T(0)} + l (b ' ij 



Proof. The first part is a consequence of Lemma 6.2 and stationary arguments. 
The second part is a consequence of Lemma 4.2, since 

the last limit being true by Lemma 6.2. □ 

Remark. Lemma 6.1 allows to extend the previous result to the cases where at 
least one of the two functions log(l ± X(p,0)) is integrable. If, for example, log(l + 
X(p,0)) e L x (0, l)andlog(l-X(p,0)) g L^O, 1), ifwetakeX(")(p,0) =min{X(p,0), 
1 - 1/n}, then X(")(p, t) < X(p, t), Vt > 0, Vp e (0, 1). Monotone Convergence The- 
orem states lim„T(™)(0) = +oo, i.e. lim„m^ = 1- Thus, by monotonicity argument 
of Lemma 6.1, moo = 1- 

With this remark in mind, we show now a "counterintuitive" example. Let 



fo(x) 



if-l<i<-l + e 
,T^( i- £l og(i-,) ) 2 if0<x<l 
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and hence 



F (x) 





±=*{l + x) 
1 -e 

1 — e log(l— x) 

1 



which corresponds to 



X (P) 



if x < -1 

if - 1 < x < -1 + e 
if - 1 + e < x < 
if < x < 1 
if 1< x 



if < p < 1 - e 



l-exp(- T ^ + i) ifl-e<p<l 



With this data, log(l + X(p, 0)) G L^O, 1) but log(l - 0)) £ I-^O, 1); 1 - e initial 
mass is close to —1 while the asymptotic solution is a 5\. 

7. Numerical examples. The analysis of the previous section left open the 
problem of the identification of the steady state in the case of the magnetization. 
Here results can be achieved only by numerical simulation of the spreading process. 
To test the numerical method, we will first derive the (explicit) solution to the pure 
drift linear equation of spreading considered in [11] as the mean field limit of the 
Sznajd model [14]. This equation reads 



9 f 9 ( ft 2w\ 



(7.1) 



namely equation (1.1) without the presence of the mean m(t). In terms of the quantile 
function X(p,t) equation (7.1) takes the form 



dX(p,t) 
dt 



- 7 X(p,t)(l-X 2 (p,t)). 



(7.2) 



Let us set 7 = —1 (spreading), and let X n (p) denote the initial datum. For any 
given p e [0, 1], equation (7.2) is an ordinary differential equation which can be easily 
integrated to give 



X(p,t) 



Xl{p)+Xl(p)e^) 1/2 ' 



(7.3) 



The asymptotic behavior of equation (7.2) can be easily deduced from the explicit 
solution. In fact, the solution converges exponentially in time to —1 if Xo(p) < 0, 
while it converges to +1 if X (p) > 0. Solution (7.3) can be inverted by using the 
definition of X(p, t). Let F (x) x e J be the initial distribution function, then, since 
X (F (x))=x 



X(F (x),t) 



(l-x 2 +x 2 e 2t ) 1/2 ' 



(7.4) 



Thus, since the function on the right of (7.4) is increasing with respect to the variable 
x, we can invert it to obtain 



X Fn 



y 



( (1 _ y 2 )e 2 t+y 2)l/2^ 



y- 



(7.5) 
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Finally, equation (7.5) finally implies 



F(y, t)=F \ y - ^ . (7.6) 

Differentiating with respect to y we conclude that, if fo(x) x S X is an initial density 
for equation (7.1), the solution in time is given by 

f(x,t) = — m f ( — ] . (7.7) 

((1 - x 2 )e 2t + x 2 f' 2 \ ((1 - x 2 )e 2t + x 2 ) 1/2 / 



- concentration T = 2 
-spreading T = 2 





Fig. 7.1. Benchmark case: evolution of density function (left hand side figure) and comparison 
between analytical and numerical solution for the quantile function (right hand side figure) 

The behavior of (7.7) shows the formation of two peaks in correspondence to the 
extremal points ±1, while in all other points of the interval X there is exponential 
decay to zero. We remark that equation (7.1) has already been solved in [11], even if 
the solution proposed in that paper is clearly wrong, as one can easily argue from the 
fact that it does not satisfy the mass conservation property. 

Using the same procedure as above, we can easily solve the problem in the opposite 
case of concentration, where 7=1. In this case, if Xq(p) denote the initial datum, 



X(p,t) 



X (p)e- 



(l-Xi(p)+X 2 (p)e-^) 1/2 ' 



(7.8) 



The solution now converges exponentially in time to zero, except for p values for which 
Xo{p) — ±1, where it remains constant. In the original formulation, the solution 
f{x,t) corresponding to the initial density fo{x) x G X is 



f(x,t) = 



-21 



((1 -x 2 )e~ 2t +x 2 ) 



,2^/2 



/o 



((1 - x 2 )e~ 2t + x 2 ) 



(7.9) 
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Note that, except for x = 0, f(x,t) converges exponentially to zero. If /o(0) > 0, the 
solution shows the formation of a peak in x = 0. 

We perform numerical simulation for different initial data in the general nonlinear 
case. First, we assume an initial symmetric datum as a benchmark (see Figure 7.1), 
where 




Cq(1 -a; 2 )(0.64-x 2 ) 1 - 3 if \x\ < 0.8; 



In this case we have m(t) = for all time t. Then we know the exact solution in order 
to perform comparison with numerical results. In Figure 7.1 we show the behavior 
of exact solution and numerical one. The last one is obtained by using standard stiff 
Runge— Kutta methods which is justified by our theoretical constructive result stated 
in the proof of Theorem 3.3. As one can see in Figure 7.1 we have a good agreement 
between analytical and numerical solution. 




■3 




Fig. 7.2. Plot of the behavior of numerical solution with symmetric initial data for concentration 
case (left hand side figure), and for spreading case (right hand side figure) 



In Figure 7.2 we sketch the plot of the quantile function X(p, t) for different times 
t both in concentration and spreading case with the same initial symmetrical data. 
As expected, in the concentration case, the limit value of all quantilcs numerically 
converge to moo = while in the spreading case the quantiles converge to —1 if 
p < .5, to 1 if p > .5. 

In Figure 7.3 we show an asymmetric case. Now, in the concentration case ^ 
given in (6.1) and quantiles converge to it. In the spreading case we can numerically 
estimate the value p^i, see (5.1), for which p < p_ x implies X(p,t) — > —1, p > p_ x 
implies X(p,t) — * 1. 

Finally, in Figure 7.4 we show a numerical "controintuitive example" analogous 
to the example given at the end of the previous section. Here, the initial datum 
consists of (1 — 2e) mass very close to —1, e mass close to 0, and the remaining e mass 
concentrated in +1. In the concentration case, 7 = 1, the mass goes asymptotically 



FIRST-ORDER CONTINUOUS MODELS OF OPINION FORMATION 



15 



■3 




■3 




Fig. 7.3. Plot of the behavior of numerical solution with non symmetric initial data for con- 
centration case (left hand side figure), and for spreading case (right hand side figure) 




---*(>, 30) 



■3 




Fig. 7.4. Plot of the behavior of a numerical "controintuitive" solution: concentration case 
tends to 8i (left hand side figure), while 1 — 2e mass goes to —1 in the spreading case (right hand 
side figure) 



to +1 while in the spreading case, 7 = —1, only 2e goes asymptotically to +1. We 
point out that the spreading case is really controintuitive. In fact, if we take X(p,T), 
with T7I, solution of the concentration case as initial data for the spreading one, the 
dynamic is as follows. Initially, (1 — e) mass decays, then it splits into two parts: 
the right one goes to +1, the left one goes to —1. The initial data can be chosen 
as close as we want to the distribution S±. Then, starting from this data (for which 
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T{t) = +00, sec (6.1)), a (1 — 2e) mass will reach —1. We note that the asymptotic 
state for spreading phenomena seems unpredictable for such concentrated initial data. 

Remark. The method we used to solve equation (7.1) represents in various case a 
possible alternative to more known methods (like the method of characteristics) able 
to reckon the solution to one-dimensional first-order partial differential equations of 
the form 

Our analysis is possible in all cases where the ordinary differential equation 




= ~<t>(X) 



is explicitly solvable. 

8. Conclusions. We investigated in this paper the spreading and/or the concen- 
tration of opinion in an organized society by means of a first-order nonlinear partial 
differential equation recently introduced in [16]. The presence of the nonlinearity ren- 
der it difficult to treat analytically the spreading case, and suitable numerical methods 
are discussed, able to capture the large-time behavior of the solution in this case. This 
work represents a first attempt for continuous approach to the formation of opinion 
in a community of agents. More complete models can be obtained by considering in 
addition the (linear or nonlinear) diffusion, which allows for a continuous steady state 
distribution function. Related problems in presence of diffusion are presently under 
study. 



(7.10) 
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